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Hadronic Vacuum Polarization Contribution to 
g — 2 of the Leptons and a(Mz) 

F. Jegerlehner a * 

a DESY-IfH Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany 

We review and compare recent calculations of hadronic vacuum polarization effects. In particular, we consider 
the anomalous magnetic moments g — 2 of the leptons and a(Mz) , the effective fine structure constant at the 
Z-resonance. 



1. VACUUM POLARIZATION AND 
CHARGE SCREENING 

Typically, charged particles in a collision of im- 
pact energy E interact electromagnetically with 
an effective charge which is the charge contained 
in the sphere of radius r ~ 1/E around the par- 
ticles. As illustrated in Fig. |l| for one of the par- 
ticles, the effective charge, due to vacuum polar- 
ization by virtual pair-creation, is larger than the 
classical charge which is seen in a large sphere 
(r — > oo). This charge screening is a particular 
kind of charge renormalization. 



m^E/2. While the lepton contributions can be 
easily calculated in QED perturbation theory the 
contribution of the strongly interacting quarks is 
not so easy to obtain. This is the issue of our 
discussion. 




7*—* e + e , fi + ft , t + t , un, dd, ■ ■ ■— > 7* 




Figure 1. 
creation. 



Vacuum polarization by virtual pair 



Commonly, Fig. ^ is represented by a Feyn- 
man diagram like the one in Fig. ^ contributing 
to muon scattering. Not surprising, the effective 
fine structure "constant" a(E) appears in many 
places in physics whenever the typical energy of 
a process is not in the classical regime. The ma- 
jor contribution to charge screening comes from 
light charged particle-antiparticle pairs of mass 



Figure 2. Feynman diagram describing the vac- 
uum polarization in muon scattering. 

1.1. Formal definition: 

The effective QED coupling constant at scale 
a/s may be written as 

(1) 



w 1 - Aa(s) 
with 

Aa(s) = -47raRc [il^s) - n^( 



(2) 



where II' (s) is the photon vacuum polarization 
function 
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, / d 4 xe^(0\Tjt; m (x)f em (0)\0) 
and j£ m (x) is the electromagnetic current. 



(3) 
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1.2. Contributions: 



with N r = 3 the color factor. 



• Leptons 

Their contribution is calculable in perturbation 
theory. The free lepton loops are affected by small 
electromagnetic corrections only. In leading order 
one obtains: 

Aoieptons(s) = 

= xTkM-/™?)- 5+° ("»?/•)] ( 4 ) 
= 0.03142 for s = Af| 



where /3e = \J\ — Amf/s. 
• Quarks 

The contribution of the five light quarks 
(it, d, s, c, b) is not calculable in perturbation the- 
ory. The free quark loops are strongly modified 
by strong interactions at low energy. The way out 
is the following: unitarity and the analyticity of 
n^(s) allow us to write (Cabbibo, Gatto 1961 1|) 

R(s') 



A 



(5) 
// mi i 



,(s) = ~V da' ; ( T ' (5) 

37T J 4m 2 S'{S' - S) 



where 



R(a) 



_ a(e + e~ 



hadrons) 



47TQ 2 (s) 

3s 



= 127rlmn' (s) 



has been measured in e + e _ -annihilation experi- 
ments up to energies above which we may calcu- 
late it in perturbative QCD. In other words 

• R{s) in known for low and moderately high 
s from the measurement of the total hadronic 
cross-section e + e~ — > hadrons. This is of par- 
ticular importance in the low energy region and 
for the resonance regions where non-perturbative 
physics comes into play. 

• R(s) may be reliably calculated for large 
s in perturbative QCD by virtue of the asymp- 
totic freedom of QCD. The leading contribution 
is given by the sum over the squares of the charges 
Q q of all quarks q 

R(s)~N c ^Q2(l + 0(a s A)) 



We thus may evaluate the hadronic part of 
the vacuum polarization by utilizing e + e _ -data 
(non-perturbative) for ^fs £ E cut ~ 40GeV plus 
the perturbative tail. Of course the data exhibit 
experimental uncertainties which will allow us to 
estimate this contribution with limited accuracy 
only. 

Before we continue to discuss the evaluation 
of the dispersion integral (||), let us remind the 
reader that the dispersion integral representation 
derives from basic properties valid for any quan- 
tum field theory: 

□ Unitarity implies the optical theorem: The 
imaginary part of the forward scattering ampli- 
tude of an elastic process A + B ^ A + B is pro- 
portional to the sum over all possible final states 
A + B -> "anything" (see Fig. g) 



Im T 



forv 



(A- 



B -> A + B) = 
y/\(s,ml,ml) a to t (A + B — > anything) 




Figure 3. Optical theorem for scattering. 



The corresponding relation for the photon 
propagator reads (see Fig. f|) 

1 

*^ S) ~ 12V 



imn;(s) = T ^ z i?(s) 




Figure 4. Optical theorem for propagation. 

□ Causality implies analyticity which may be 
expressed in form of a so-called (subtracted) dis- 
persion relation 



n;( g 2 )-n'(o) 



IT 



ds 



imn;( s ) 

s (s — q 2 — ie) 
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The latter, together with the optical theorem, di- 
rectly implies the validity of (|h. Note that its 
validity is based on general principles and holds 
beyond perturbation theory. It is the basis of all 
non-perturbative evaluations of hadronic vacuum 
polarization effects in terms of experimental data. 



1.3. Description of the 



" data. 



In a recent reanalysis Q| of the hadronic vac- 
uum polarization we have collected all published 
e + e~-data plus some unpublished data (see also 
j3|). The data sets obtained for different energy 
regions have been displayed in a number of figures 
in Ref. [Q, and we will not repeat them here. The 
following integration procedure has been used 
for the evaluation of the dispersion integral: 

O Take data as they are and apply trapezoidal 
rule (connecting data points by straight 
lines) for integration. 

© To combine results from different experi- 
ments: i) integrate data for individual ex- 
periments and combine the results, ii) com- 
bine data from different experiments before 
integration and integrate combined "inte- 
grand" . Check consistency of the two pos- 
sible procedures to estimate the reliability 
of the results. 

© Error analysis: 1) statistical errors are 
added in quadrature, 2) systematic errors 
are added linearly for different experiments, 
3) combined results are obtained by taking 
weighted averages. 4) all errors are added 
in quadrature for "independent" data sets. 
We assume this to be allowed in particular 
for different energy regions and/or different 
accelerators. 

Resonances 

have been parametrized by Breit-Wigner 
shapes with parameters taken from the Par- 
ticle Data Tables {§. 

Our recent update (Eidelman, Jegerlehner 1995) 

^ { hLons(M 2 z ) = 0-0280 ± 0.0007 

together with the leptonic term (0) yields an 



Table 1 

Comparison of 



estimates of Aaj^Mf ) 



Author (Year) [Ref.] 



0.0285 (7) Jegerlehner (86) 

0.0283(12) Lynn et al. (87) LC 

0.0287 (9) Burkhardt et al. (89) f) 

0.0282 (9) Jegerlehner (91) § 

0.0273 (4) Martin, Zeppenfeld (94) @ 

0.0280 (7) Eidelman, Jegerlehner (95) 

0.0280 (7) Burkhardt, Pietrzyk (95) fl 

0.0275 (5) Swartz (95) Jll| 

0.0289 (4) Adel,Yndurain (95) || 



effective fine structure constant at the Z- 
resonance: 

a-\M z ) = 128.89 ±0.09 . 
1.4. Results for a(M z ): 

In Table § and Fig. | we show a comparison of 
results obtained by different authors. All values 
have been rescaled to M z = 91.1887 GeV. The 
small top quark and the W contributions are not 
included in a{M z ). For a comparison of the ear- 
lier results |lj| we refer to |i| . The differences are 
mainly due to a different treatment of the system- 
atic errors and/or different model assumptions. 
Detail of our analysis are given in Table 2. Note 
that a reduction of the error would require a pre- 
cision measurement of the <7tot(e + e - — > hadrons) 
from 1 GeV to about 10 GeV. 



Jegerlehner (86) 
Lynn et al. (87) 
Burkhardt et al. (89) 
Jegerlehner (91) 
Swartz (94) 

Martin & Zeppenfeld (94) 
Eidelman & Jegerlehner (95) 
Burkhardt & Pietrzyk (95) 
Swartz (95) 
Adel & Yndurain (95) 



128.6 128.7 128.8 128.9 129.0 129.1 129.2 129,3 



Figure 5. Different estimates of a(M z ) 



-l 



A few remarks concerning the results presented 
in Tab. |l| are in order here. Updates seemed to be 
justified from time to time when new experimen- 
tal data became available. While there have been 
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Table 2 



Contributions to Aa^A x 10 4 and relative(rel) and absolute(abs) errors in percent. 



final state 


range (GeV) 


contribution 




contribution 


rel % 


abs % 


P 


(0.28, 0.81) 


26.08 


(0.29) 


(0.62) 












(0.42, 0.81) 


2.93 


(0.04) 


(0.08) 












(1.00, 1.04) 


5.08 


(0.07) 


(0.12) 


Resonances: 


46.61 (1.08) 


2.3 


0.4 






11.34 


(0.55) 


(0.61) 










T 




1.18 


(0.05) 


(0.06) 










liadrons 


(0.81, 1.40) 


13.83 


(0.15) 


(0.79) 


Background: 








hadrons 


(1.40, 3.10) 


27.62 


(0.32) 


(4.01) 


E < Mj^ 


41.45 (4.11) 


9.9 


1.5 


hadrons 


(3.10, 3.60) 


5.82 


(0.30) 


(1.12) 


Mj fil < E < 3.6 GeV 


5.82 (1.16) 


19.9 


0.4 


hadrons 


(3.60, 9.46) 


50.60 


(0.24) 


(3.33) 


3.6 GeV < E < M r 


50.60 (3.34) 


6.6 


1.2 


hadrons 


(9.46, 40.0) 


93.07 


(0.86) 


(3.39) 


Mr < E < 40 GeV 


93.07 (3.50) 


3.8 


1.2 


QCD 


(40.0, oo ) 


42.82 


(0.00) 


(0.10) 


40 GeV < E QCD 


42.82 (0.10) 


0.2 


0.0 


total 




280.37 


(1.18) 


(6.43) 


total 


280.4 (6.54) 




2.3 



very little truly new experimental results, some 
groups have published updated results which were 
available before as preliminary data only. Exam- 
ples are the ND results @ up to 1.4 GeV, the 
DM2 results jl| between 1.35 and 2.3 GeV and 
the Crystal Ball results |l6| between 5.0 and 7.4 
GeV. New results from VEPP-4 @, from 7.2 to 
9.5 GeV, have been included as well. In addition 
in |l|] we have made an effort to collect as much 
as possible all the available data. 

Additional motivations for performing the up- 
date II were the following: the recent issues of 
the Review of Particle Properties Q had some 
improvements on the resonance parameters and 
there was progress in calculating R(s) at high 
energies in perturbative QCD. For the muon 
anomaly the study of the low energy end by 
means of chiral perturbation theory allowed us 
to reduce potential model-dependences of earlier 
approaches. 

The update Q| attempted to estimate an "of- 
ficial" number which could be used by LEP ex- 
periments in their analyses. As compared to ||, 
a more conservative guess for the uncertainty in 
the continuum above the p up to the J/^> was 
made, and therefore, the uncertainty increased 
from 0.0007 to 0.0009. Later, the Crystal Ball 
(CB) collaboration carefully reanalyzed their old 
e + e~- annihilation data and obtained R(s) val- 
ues substantially lower than the Mark I data and 
in agreement with other experiments (LENA). 



The results were in much better agreement with 
perturbative QCD. The change of the data was 
mainly due to an up-to-date treatment of the 
QED radiative corrections and r-subtraction. In 
Ref. H I included the updated results from CB 
and discarded the Mark I data, which system- 
atically lie 28% higher, on average. Note that, 
because of the much larger systematic uncertain- 
ties of the Mark I data, the results are affected in 
a minor way if we exclude (50.60(0.24) (3.33)) or 
include (50.79(0.20)(3.20)) the Mark I data. The 
inclusion of the CB data lead to a reduction of 
previous results by -0.0005. While in Refs. [|J^J| 
the p region was parametrized by an improved 
Gounaris-Sakurai formula, in our last update |]J 
we apply trapezoidal rule integration and include 
the final data from VEPP-2M OLYA, ND and 
CMD up to 1.4 GeV. This extended analysis con- 
firmed our previous results. From 1.35 to 2.3 GeV 
new data from ORSAY DCI DM2 lead to a re- 
duction of the error in this region. The error now 
agrees with the earlier estimate in ||. 

While the analysis jl(| is similar in spirit to 
ours and reproduced our result, the estimate || is 
biased by believing in perturbative QCD down to 
3 GeV; some experimental data used in the J /if? 
and T resonance regions are rescaled using per- 
turbative QCD. In contrast to this estimate which 
yields a low value, the estimate |l2] ], which also 
relies on perturbative QCD, finds a large value; 
both estimates result in much lower errors than 
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what can be justified by the data alone. In the 
analysis of data are fitted before integration. 
This requires a guess of the functional form of the 
integrand which is then fitted to the data. In ad- 
dition one has to assume some kind of correlation 
matrix between the data points. One consequence 
of the method applied is that the inclusion of one 
additional unpublished data point |l6| ] led to a 
shift of the result by 1 a. 

In Fig. U we show the result for a{— Q 2 ) as a 
function of E = \Q\ for low t = —Q 2 in the space- 
like region, which is relevant for the ^-channel 
contribution to Bhabha scattering. Note the 



0.06- 

0.05 

0.04 

0.03 

0.02 

0.01 

0.00 



leptons + 5 quarks 
leptons only 



-30.0 -25.0 -20.0 -15.0 -10.0 -5.0 0.0 
E (GeV) 

Figure 6. Aa(- Q 2 ) in the spacelike region. 
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Figure 7. Uncertainty of Aa(t) 



of a(— Q 2 ) is possible at LEP by an appropri- 
ate analysis of the small angle Bhabha scattering 
data. The values of Q 2 should be low enough such 
that the i-channel contribution is clearly domi- 
nant. In Fig. ^ we show the uncertainty in % of 
our evaluation of the running charge in the low 
energy region. 



1.5. Parametrizations. 



We have performed fits of Aa{^ d (s). The 
parametrizations^ and the best fit parameters we 
found are the following: 

For -(2GeV) 2 < s < (0.25GeV) 2 the best fit is 



Ci < In II - c 2 s| + c 2 - 



c 3 s 



h 1- c 4 (s/so) 

c 3 - s 



C3 - s 

2 



with s = -(2GeV) 2 , h = 9.3055 x 1(T 3 GeV~ 2 , 
ci = 2.2694240 x 10~ 3 , c 2 = 8.073429GeV" 2 , 
c 3 = 0.1636393GeV 2 , c 4 = -3.35455 x 10~ 5 . 

In the range -(20GeV) 2 < s < -(2GeV) 2 we find 



A4 5 a } d (s) 



Ci < In II - c 2 s\ + c 2 - 



c 3 s 



C3 



-h 1- c 4 (s/so) 

c 3 - s 



with s = -(20GeV) 2 , h = 9.3055x lO" 3 GeV" 2 , 
ci = 2.8668314 x 10" 3 , c 2 = 0.3514608GeV" 2 , 
c 3 = 0.5496359GeV 2 , c 4 = 1.989233 x 10~ 4 . 

At higher energies, in the ranges E m - m < E < 

Em&x, we have 



Ci + c 2 (ln(s/s ) 

+c 3 (so/s- 1) + c 4 ((s /s) 2 



1)) 



dramatic increase of the effective charge at low 
spacelike momenta. This shows that the classi- 
cal limit is difficult to reach in a scattering ex- 
periment. A clean measurement of the running 



2 Fortran routines hadr5.f and alphaQED.f for calcu- 
lating Att [ 5 ^(s) and a(s), respectiv ely, are available from 



WWW http://www.ifh.de/ fjeger/ 
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where s = E\E\ and s = E \E \ with 



Exnin 


250 


80 


40 




1000 


250 


80 


£■0 


1000 


91.1888 


80 


ci x 10^ 


4.2092260 


2.8039809 


2.7266588 


c 2 x 10 3 


2.9233438 


2.9373798 


2.9285045 


c 3 x 10 3 


-0.3296691 


-2.8432352 


-4.7720564 


c 4 x 10 4 


0.0034324 


-5.2537734 


7.7295507 


-^min 


-200 


-1000 






-20 


-200 






-100 


-1000 




ci x 10^ 


2.8526291 


4.2069394 




c 2 x 10 3 


2.9520725 


2.9253566 




c 3 x 10 3 


-2.7906310 


-0.6778245 




c 4 x 10 4 


0.6417453 


0.0932141 





At the higher energies, E > 40 GeV, we 
may fit the hadronic contribution using the fol- 
lowing effective parameters in a 0{ct s ) pertur- 
bative QCD formula: effective quark masses: 
m u<d , s , c<b = 0.067,0.089,0.231,1.299,4.500 GeV 
together with an effective a s = 0.102. Thus 

(h(y d ) + h(y s ) + h(y b ) + 4(% u ) + h(y c ))) 



where 
h(y) = 

g{y) = 

with yi 



5/3 + »-(l + y/2) «/(») 



4m? /s. 



1 arctan(l/ v / y — 1) 



for 
for 



y > i 
y < i 



which may be evaluated in terms of the experi- 
mentally accessible quantity 



R(s) = CTtot , ( + e ' 



*hadrons) 




hadrons 

Figure 8. Leading vacuum polarization contribu- 
tion to g — 2. 

The kernel K(s) of (§), depicted in Fig. [|, is a 
smooth bounded function. 



K(s) 
1 

0.63.. 




Figure 9. The kernel K(s) as a function of energy. 



2. VACUUM POLARIZATION CONTRI- 
BUTION TO 5 2 

The leading hadronic contribution to g — 2 is 
given by the diagram Fig. ||, where the blob repre- 
sents the one-particle irreducible contribution to 
the photon propagator. To evaluate the diagram 
Fig. H we can make use of the representation (||) 
for the vacuum polarization function. Again one 
obtains a dispersion integral (Bouchiat, Michel 
1961, Durand III 1962, Gourdin, de Rafael 1969) 



4m? 



Note that the energy denominator 1/s 2 — 1/E 4 
dramatically enhances the non-perturbative low 
energy contribution 

2.1. Results for a*f d : 

The present status od the hadronic vacuum 
contribution to g — 2 is summarized in Table || 
and Fig. [l^. Only the leading hadronic vacuum 
polarization diagram is included. Our estimate is 

fl had = 702.35(15.28) x 10~ 10 . 

Our values for a^ ad and a^ ad are given in Table 
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Table 3 



a had . ^giu anc j re lative(rel) and absolute(abs) errors in percent. 



final state 


range (GeV) 




a had . 10 io 


rel % 


abs % 


P 


(0.28, 0.81) 


426.66 ( 


5.61) 


10.62) 


2.8% 


1.7% 




(0.42, 0.81) 


37.76 


( 0.45) 


( 1-02) 


3.0% 


0.2% 


<i> 


(1.00, 1.04) 


38.55 


( 0.54) 


( 0.89) 


2.7% 


0.1% 


J/i> 




8.60 


[ 0.41) 


( 0.40) 


6.7% 


0.1% 


T 




0.10 


[ 0.00) 


( 0.01) 


6.7% 


0.0% 


hadrons 


(0.81, 1.40) 


112.85 


( 1.33) 


( 5.49) 


5.0% 


0.8% 


hadrons 


(1.40, 3.10) 


56.43 


( 0.45) 


( 7.22) 


12.8% 


1.0% 


hadrons 


(3.10, 3.60) 


4.47 


( 0.23) 


( 0.86) 


19.9% 


0.1% 


hadrons 


(3.60, 9.46) 


14.06 


( 0.07) 


( 0.90) 


6.5% 


0.1% 


hadrons 


(9.46, 40.0) 


2.70 


' 0.03) 


( 0.13) 


4.9% 


0.0% 


QCD 


(40.0, oo ) 


0.16 


{ 0.00) 


( 0.00) 


0.2% 


0.0% 


total 




702.35 ( 


5.85) 


14.09) 


2.2 % 


2.2% 



Table 4 

Various estimates of a'' ad x 10 10 



a'f d x 10 



10 



Author (Year) [Ref.] 



650 (50) 
680 (90) 
663 (85) 
702 (80) 
684 (11) 
707 (6)(16) 
710(10)(5) 
705 (6) (5) 
724 (7)(26) 
699 (4) (2) 
702 (6)(14) 
711 (10) 
702 (8) (13) 



Gourdin, de Rafael (69) |2 
Bramon, Etim, Greco (72) 
Barger, Long, Olsen 
Calmet et al. (78) 
Barkov et al. (85) |2f 
Kinoshita et al. (85) 
Casas et al. (85) pq ] 
Martinovic, Dubnicka (90) |2|] 
Jegerlehner (91) ]30[ | 
Dubnickova et al. (92) |l| 
Eidelman, Jegerlehner (95) 
Adel, Yndurain (95) HI 




Worstell, Brown (95) 32] 



— Bouchiat k Michel (1961) 
-b Kinoshita & Oakes (1967) 

Gourdin & de Rafael (1969) 

— Bramon, Etim k Greco (1972) 
h Barger, Long k Olsen (1975) 

— Calmet et al. (1976) 

— Narison (1978) 
Barkov et al. (1985) 
Kinoshita et al. (1985) 
Casas et al. (1985) 
Martinovic k Dubnicka (1990) 

i Jegerlehner (1991) 
Dubnickova et al. (1992) 
Eidelman k Jegerlehner (1995) 
Adel k Yndurain (1995) 
Worstell k Brown (1995) 



1 10.00 6(0.00 680.00 700.00 720.00 740.0 



Figure 10. Evaluations of aj^ 
ferent authors. 



hadr 



10 10 by dif- 



3. CONCLUSION AND OUTLOOK 

Muon anomalous magnetic moment: The 
present uncertainty due to hadronic vacuum po- 
larization is 



(.9 ~ 2) A 



5a u 



16 x 10 



-10 



This can be reduced to 8 x 10 by ongoing mea- 



surements at VEPP 2M in Novosibirsk |35| and 
in a forthcoming experiment at DA<I>NE in Fras- 
cati pa]. This has to be compared with the ex- 



perimental precision of the forthcoming BNL ex- 
periment E821 at Brookhaven j33|. The present 
experimental precision is 84 x 10~ 10 was obtained 
in the CERN experiment @. A look at Table | 



shows that a measurement of e + e 



hadrons 



for 1.5 GeV to 3.6 GeV to 1% accuracy is ur- 
gently needed in order to be able to fully exploit 
the BNL measurement of a^. 
""»■ Bhabha scattering: the hadronic uncertainty 
is not a limiting factor in small angle Bhabha 
scattering (luminosity monitoring). The esti- 
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Table 5 

Contributions to a|? ad • 10 14 and a^ ad ■ 10 8 . 



final state 


range (GeV) 


of d • 10 i4 


(stat) (syst) 


af d • 10 8 


(stat) (syst) 


P 


(0.28, 0.81) 


118.57 ( 


1.59) ( 2.97) 


137.42 ( 


1.61) ( 3.33) 




(0.42, 0.81) 


10.07 ( 


0.12) ( 0.27) 


14.62 ( 


0.18) ( 0.39) 




(1.00, 1.04) 


9.86 ( 


0.14) ( 0.23) 


20.36 ( 


0.29) ( 0.47) 


J/i/> 




2.04 ( 


0.10) ( 0.09) 


12.80 ( 


0.60) ( 0.61) 


T 




0.02 ( 


0.00) ( 0.00) 


0.24 ( 


0.01) ( 0.01) 


hadrons 


(0.81, 1.40) 


29.16 ( 


0.35) ( 1.41) 


56.22 ( 


0.62) ( 3.00) 


hadrons 


(1.40, 3.10) 


13.70 ( 


0.11) ( 1.75) 


56.32 ( 


0.54) ( 7.60) 


hadrons 


(3.10, 3.60) 


1.06 ( 


0.06) ( 0.20) 


6.66 ( 


0.35) ( 1.28) 


hadrons 


(3.60, 9.46) 


3.31 ( 


0.02) ( 0.21) 


26.43 ( 


0.13) ( 1.70) 


hadrons 


(9.46, 40.0) 


0.63 ( 


0.01) ( 0.03) 


6.78 ( 


0.06) ( 0.33) 


QCD 


(40.0, oo) 


0.04 ( 


0.00) ( 0.00) 


0.45 ( 


0.00) ( 0.00) 


total 




188.47 ( 


1.65) ( 3.75) 


338.30 ( 


1.97) ( 9.12) 



mated uncertainty is ~ 0.06 % . 
"<+■ Precision measurements at LEP and SCL: the 
hadronic uncertainty is close to become the lim- 
iting factor for precision measurements of some 
observables like 



sin 2 6 



c rr 



(5sin 2 e ~ 0.00025 



This compares to the present accuracy of 
S sin 2 6 - 0.00030 of LEP and SLC data. 

In principle one could try to do a direct mea- 
surement of (g — 2) M h a( j and a(Mz)- Unfortu- 
nately, in this case one would lose the most impor- 
tant predictions and in addition the result would 
be model dependent !. 

Reliable theoretical calculations are not in 
sight. 

A real breakthrough in improving on the 
hadronic uncertainty will require dedicated efforts 
in high precision measurements of R(s) in a wide 
energy range. 

4. APPENDIX: Analytic expressions. 

In many cases, in some interval 4m 2 < s% < 
S2 < oo, one may approximate R(s) by simple 
analytic expressions and the dispersion integral 



Aa had {s) 



At: 2 a 



V ds 



I <7had{s') 

s' — s 

1 1 



= f^i <-- 



R(s') 



may be performed analytically. 

We consider a few examples in the following. 

• Linear extrapolation 

For example, for non-resonant contributions 
(neighboring) data points may be approximated 
by a linear function 

R{s) = a + b^s 

in the cm. energy which gives a contribution^] 

Aa had (s) = _|2aln— -(a + ^ln^-^ 

3tt I Wi si — s JK 'V 



where Wi = y/sl, W = y/s and 



f(m = i ^N^I-H™l);M^ 2 >o 

I(W) \2|W|(arctanf" 



arctanJ^P); q 2 = W 2 < 0. 



Wi 



• Narrow width resonance 



The contribution from a zero width resonance 

127T 2 

cr N w(s) = -r^— T e+e -S(s - M R ) 

JMR 

3 Eq. (7) of Ref. H contains several misprints 



9 



or 



Rnw{s) = ^T e+e -S(s- M%) 



is given by 
Aa res (s) = 



which in the limit \s\ ^S> Ml becomes 

A, , 3r„ 

V ; ttMji 
• Breit-Wigner resonance 

The contribution from a Breit-Wigner resonance 

OBw{s) ■ 

or 



ill 

4 



3tt rr e+e ~ 

a (^-M R ) 2 + 

9 rr e+e - 

4a 2 (^_ Mr )2 + ^ 



47ra 



{/(O)-J(W)} 



is given by 
Aa res (s) = 

where 
I(W) = i 

r 2lc / \ 

\ VK-Mr— tc \ Ln Wi-W 111 VCi-AfR-to y 

1 A„ W2+W _ t W 2 -M n -ic ' 

W+M R +ic \^ Wi+W Wt-MR-ic 

-/l.C.} 

with c = r/2. This may be written as 



Aa res (s) 
{s(s 
1 7r — arctan 



3r„ 



My. 



M 2 R + 3c 2 



■kcxMr M 2 +c 2 (s-M 2 +c 2 ) 2 + M 2 r 2 

arctan ■ 



W 2 -M R 



Mft—Wi 



r 

4M n 



3(s _ 3 A f 2 +c 2 )ln (^ = M.) a +c 



si — s I 



+ (a + M 2 + C 2 )(M 2 + c 2 ) In 
+2M fl (M| + c 2 )f(W) 
-K2(( S -M 2 +c 2 ) 2 + M 2 r 2 )ln^ 

with /(W) given above. For W x < Mr < W 2 
and r <C M/j this may be approximated by 

_ 3r e+e - g(a - M I + 3c 2 ) 
reslSj " aM R (s - Ml + c 2 ) 2 + M 2 r 2 

which agrees with Eqs. (3) and (4) in the limits 
r 2 < \s - M%\,M% and \s\ » M|, respectively. 



• Breit-Wigner resonance: Field theory 
version 

Finally, we consider a field theoretic form of a 
Breit-Wigner resonance obtained by the Dyson 
summation of a massive spin 1 transversal part 
of the propagator in the approximation that the 
imaginary part of the self-energy yields the width 
by Imlly(My) = MyTv near resonance. 



(Jbw(s) 



12tt T„ 



sT 2 



Ml 



r (s - Ml) 2 + M 2 R Y 2 



which yields 



lAUresyh) - 7rctMj? (s- M 2 R ) 2 +M 2 T 2 



{( 



7r — arctan 



si-Mi 



arctan 



TMr 
"72" — 



Mr (s-M^-r 2 ) 



f In I «a— « I In I »2--M|-iM B r A 1 



and reduces to 



Aa res (s) 



3r e+e ^ 8 ( g - M 2 - r 2 ) 
aM R (s-M 2 ) 2 + M 2 r 2 



for si <C Ml <C S2 and T <C Mr. Again we have 
the known limits for small V and for large \s\. 
For broad resonances the different parametriza- 
tions of the resonance in general yield very dif- 
ferent results. Therefore, it is important to know 
how a resonance was parametrized to get the res- 
onance parameters like Mr and T. For narrow 
resonances as considered in our application re- 
sults are not affected in a relevant way by using 
different parametrizations. 
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